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One of the most challenging issues in particle physics is to study QCD in extreme conditions. Precise de- 
termination of the QCD phase diagram on temperature T and chemical potential fi plane will provide valuable 
information for quark-gluon plasma (QGP) and neutron star physics. We present results for phase structure on the 
(/i, T) plane for lattice QCD with Wilson fermions from strong coupling Hamiltonian analysis and Kogut-Susskind 
Fermions from Lagrangian Monte Carlo simulations at intermediate coupling. 



1. Introduction 

For two massless quarks, several QCD inspired 
models suggest the existence of a tricritical point 
on the (jx, T) plane separating the first order tran- 
sition line at lower T and larger \i, and the sec- 
ond order transition line at higher T and smaller 
fi. Lattice gauge theory (LGT) should in prin- 
ciple give more reliable information on the QCD 
phase diagram. Heavy-ion collision experiments 
may determine whether the theoretical prediction 
is correct. 

In Lagrangian formulation of SU(3) LGT, the 
chemical potential fj, is introduced by making 
the replacement: 1/4(2;) = e^ a U^{x), U±(x) — 
e~^ a U^(x). It is well known that, for real fi, 
the effective fermionic action is complex and the 
standard MC technique does not apply. The re- 
cent years have seen enormous efforts £Q on solv- 
ing this problem, and some very interesting in- 
formation |2I3I4| on the phase diagram for QCD 
with KS fermions at large T and small fx has been 
obtained from MC simulations. There is also a 
strong coupling analysis of the color supercon- 
ductivity phase at T = 0. 

Nevertheless, to precisely locate the tricritical 
point and critical line at large \i is still an ex- 
tremely difficult task. QCD at large \i is of partic- 
ular importance for neutron star physics. Hamil- 
tonian formulation of LGT doesn't encounter the 
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notorious "complex action problem" . Recently, 
we proposed a Hamiltonian approach |6l7| to LGT 
with naive fermions at finite //, and extended it 
to Wilson fermions [S]. The chiral phase transi- 
tion at T = and some finite fic was found to 
be of first order. In Ref. [H], the authors studied 
lattice QCD with SLAC fermions at T = in the 
strong coupling limit. 

In this paper, we study the phase diagram on 
the (/x, T) plane by Lagrangian MC simulations 
with Kogut-Susskind quarks and Hamiltonian 
strong coupling analysis with Wilson fermions. 

2. MC Simulation with Nf = 2 and 4 

For small /z, some approximation can be used in 
simulations. If \i is purely imaginary, the effective 
fermionic action is real and the traditional MC 
methods still applies. Using Taylor expansion, 
the physical observables at imaginary \x can be 
extrapolated to real /1. 

For imaginary [i — ifii, Roberge and Weiss fTUj 
find that the partition function is periodic in 
Hi/T with periodicity 27r/3. The Polyakov loop is 
singular at fii/T = 27r(fc + l/2)/3, with k an inte- 
ger. In pure gauge theory, the partition function 
has the universal Z(3) symmetry. In full QCD, 
the fermionic determinant breaks the Z(3) sym- 
metry. When fij is increased from 0, there are 
phase transitions at fXi/T = 2n(k + l/2)/3. For 
[iifT < 7r/3, the simulation is physical and for 
Hi/T > 7r/3, the result is only a simple copy 
of that for /xj/T < 7r/3, which is the artifact 
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of imaginary /z. So this method works only for 
Mj/T < tt/3. 

We have simulated QCD with KS fermions us- 
ing the standard R and HMC algorithms, for 
Nf — 2 and Nf — 4 respectively. We modified 
the MILC code to imaginary ifii. The molecu- 
lar dynamics for each configuration consists of 20 
steps with step size St is 0.02. The lattice size 
is 8 3 x 4. The total trajectories are 4000. The 
measurements were done for every 40 trajectories. 
Fig. n an( i Fig. [5] show our results. The data for 
the chiral condensate, Polyakov loop and the its 
phase are consistent with Refs. |2"I3 . In addition, 
we measured the plaquette and fermionic energy 
density. 

3. Tricritical Point at Strong Coupling 

In the Hamiltonian formulation, one can study 
QCD at real fi. At strong coupling, the effective 
Hamiltonian with four fermion interactions |8lll] 
is the result from integrating out the gauge fields. 
For convenience, we rescale the chemical potential 
and temperature as // = fi/(3K/a) and T 1 = 
T / (3K / a) , with K being the effective coupling of 
four fermion interactions. At the second order 
phase transition where the chiral condensate and 
the dynamical mass of quark vanish continuously, 
there is only one global minimum in the grand 
thermodynamic potential. For Wilson fermions 
and Nf/N c < 1, where N c = 3 is the number of 
colors, we obtain an equation for the critical line 
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which is depicted by the dotted line for r = 1 in 
Fig. El 

Below some finite Tg, there is a first order chiral 
phase transition line 

^' c = l + 2r 2 . (2) 

from some finite Tg down to T' = 0. this is illus- 
trated by the solid line for r = 1 in Fig. [21 

The point when two lines described by Eq. (JIJ 
and Eq. @ join at lower T' is the tricritical point 



(/4 T 3 '). Forr = I, we find {^,r 3 ) = (3,0.4498), 
i.e. the circle in Fig. The phase structure for 
any r ^ is qualitatively the same. Details can 
be found in Ref. [T2]. 



4. Summary 

In the preceding sections, we present some pre- 
liminary MC data for QCD with KS fermions at 
some finite /3 and imaginary fi. These can be 
converted to some T and real fi. The results for 
Nf = 2 might indicate the existence of a second 
order chiral phase transition, if the data are ex- 
trapolated to the chiral limit. Those for Nf = 4 
indicate the existence of a first order phase tran- 
sition. To determine whether a tricritical point 
exists for Nt = 2, further detailed study is re- 
quired. 

We also investigate the QCD phase diagram 
in Hamiltonian lattice formulation with Wilson 
fermions. At the strong coupling, we find a tri- 
critical point on the (fj,, T) plane, which has not 
been found in previous work in the Hamiltonian 
formalism with KS or naive fermions. Our find- 
ings imply that on the (/x, T) plane, at least in the 
Hamiltonian formulation at the strong coupling, 
the phase structure of QCD with Wilson fermions 
(without species doubling) might be qualitatively 
different from naive or KS fermions (with species 
doubling) . 
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Figure 1. Results for KS quarks at j3 = 5.1, am — 
0.05, and Nf = 4. <fi is the phase of Polyakov loop. 
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Figure 2. Results for KS quarks at (3 = 5.3, am 
0.25, and JV> = 2. 
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Figure 3. Phase diagram for Wilson quarks. The 
solid and dotted lines stand respectively for the 
1st and 2nd order transitions. The circle is the 
tricritical point. 



